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ABSTRACT 

We construct an evolutionary spectral energy distribution (SED) model of a starburst 
region, from the ultraviolet to submillimetre wavelengths. This model allows us to 
derive the star formation rate, optical depth by dust and apparent effective radius of 
starburst regions at various wavelengths; as a result, the intrinsic surface brightness 
of starburst regions can be derived. Using this SED model, we analyse 16 UV-selected 
starburst galaxies and 10 ultraluminous infrared galaxies. The derived star formation 
rates and optical depths are compared with emission line measurements and found to 
be consistent. The derived apparent effective radii are also consistent with observa- 
tions. From the SED analysis, we find a bimodal property of the star formation rate 
with the optical depth and the compactness of stellar distributions. While mild star- 
bursts have a hmiting intrinsic surface brightness Lf,oir~'^ ~ 10^^ L0kpc~^, intense 
starbursts tend to be more heavily obscured and concentrated within a characteristic 
scale of Te ~ 0.3 kpc. We suggest that the mild starbursts can be triggered by a self- 
gravitating disc instability in which feedback is effective in the shallow gravitational 
potential. On the other hand, the intense starbursts can be induced via an external 
dynamical perturbation like galaxy merging, in which feedback is less effective due 
to the deep gravitational potential attained by the large gas concentration within the 
central starburst region. 

Key words: galaxies: starburst - dust, extinction - infrared: galaxies - ultraviolet: 
galaxies - submillimetre. 



1 INTRODUCTION 

Starbursts are an intensive mode of star formation in galax- 
ies (Storchi-Bergmann, Calzetti & Kinney 1994; Kenni- 
cutt 1998). Ultraluminous infrared galaxies (ULIRGs) with 
Lm > lO^'^ L© are probably the most active and luminous 
starburst galaxies. Indeed recent ISO observations suggest 
that many ULIRGs are powered by starbursts with little 
contribution from AGN (Genzel et al. 1998; Lutz et al. 
1998). On the other hand, starburst galaxies of lower ac- 
tivity are HII galaxies and blue compact dwarf galaxies, 
known as UV-selected starburst galaxies (UVSBGs; Kin- 
ney et al. 1993; Storch-Bergmann et al. 1994; Calzetti, Kin- 
ney & Storchi-Bergmann 1994; McQuade, Calzetti & Kin- 
ney 1995; Gordon, Calzetti & Witt 1997; Takagi, Arimoto 
& Vansevicius 1999). 



Recently, an upper limit of 'bolometric surface bright- 
ness' was given by Meurer et al. (1997) for a sample of star- 
burst galaxies observed in the rest frame 1) UV, 2) far in- 
frared (FIR) and Ha, and 3) 21 cm radio continuum emis- 
sion. This limit seems to be physically associated with an 
instability in gaseous discs (e.g. Toomre 1964; Quirk 1972) 
as Kennicutt (1989) had already suggested for the star for- 
mation in normal disc galaxies. In their study, the effec- 
tive radii are determined from observations at the different 
wavelengths for each galaxy. However, if a galaxy is opti- 
cally thick the observed effective radius strongly depends on 
the degree of extinction, since the observed light comes out 
solely from the region in which the optical depth is less than 
unity. Moreover, many ULIRGs have multiple starburst re- 
gions (e.g. Scoville et al. 2000; Surace & Sanders 2000), while 
most of UVSBGs show only one major central starburst re- 
gion. These effects suggest that the panchromatic intensity 
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limit on starbursts proposed by Meurer et al. (1997) may 
not be appropriate for ULIRGs. 

To determine the intensity limit of staxbursts with var- 
ious levels of activity ranging from UVSBGs to ULIRGs, it 
is necessary to estimate both star formation rates (SFRs) 
and effective radii in a unified manner, whatever the de- 
gree of extinction. While typical Y-band optical depths in 
UVSBGs are tv ^ 0.3 — 2 (Takagi, Arimoto, & Vansevicius 
1999; Storchi-Bergmann et al. 1994), those in ULIRGs reach 
TV ^ 5 — 50 (Genzel et al. 1998). Such a wide variation in tv 
causes a serious problem in deriving SFRs, since none of the 
currently used SFR indicators, such as Ha (or other Balmer 
lines), UV continuum, or FIR luminosities, are commonly 
applicable for starbursts ranging from UVSBGs to ULIRGs. 
The FIR luminosity is not a good measure of bolometric 
luminosity when starburst regions are optically thin, while 
Ha luminosity is difficult to use if the mean surface gas 
density becomes larger than 50 Mq pc~'^, corresponding to 
TV ~ 2, causing a significant extinction at Ha (Kcnnicutt 
1998). What is worse, by using Br7 and FIR luminosities, 
Kennicutt (1998) has shown that the FIR luminosities give 
systematically higher SFRs (by a factor of ~ 2) than the 
Br7S. Sullivan ct al. (2000) compared the Ha- and UV- 
derived SFRs of nearby starbursts to find that starburst 
galaxies are typically over-luminous in the UV for a given 
Ha luminosity. The effect is strongest in the less luminous 
galaxies. 

When comparing ULIRGs and UVSBGs, the bolomet- 
ric surface brightness should be derived for each individual 
starburst region in ULIRGs. The wealth of ground-based 
low-resolution imaging surveys have provided little informa- 
tion on the stars in the starburst region near the confusion 
limit < 1". Even the Hubble Space Telescope (HST) with 
its higher resolution < 0.1" fails to go deep enough into 
the starburst regions, which are obscured by dust at optical 
wavelengths. The Near-Infrared Camera and Multi-Object 
Spectrometer (NICMOS) on HST can observe the obscured 
starburst regions more directly at near infrared (NIR) wave- 
lengths where the effects of dust extinction are reduced sig- 
nificantly compared to visual wavelengths. Scoville et al. 
(2000) studied the morphologies of 24 ULIRGs with NIC- 
MOS and found that light profiles of nine of them were well 
fit by an r^^^ law, rather than an exponential profile. The 
apparent effective radii of ULIRGs tend to be more compact 
than the extent of gas observed with high-resolution imaging 
of the CO emission (Bryant & Scoville 1999). This suggests 
that the distribution of stars is more concentrated than that 
of dust and gas in the starburst regions. Thus, the dimming 
of light due to the surrounding dust can be so large that the 
real distribution of stars deviates from the apparent distri- 
bution of the light even the longest observable wavelength 
with NICMOS of 2.2 /xm. 

It is thus crucial to establish a measure of both SFRs 
and the spatial distribution of stars which can be applied 
to starbursts of any optical depth. Clearly, a new recipe to 
derive both SFRs and the stellar distributions is required to 
investigate the physical properties of starbursts comprehen- 
sively. 

In this paper, we will use spectral energy distributions 
(SEDs) from the far ultraviolet (FUV) to submillimetre 
(submm) to derive SFRs of starburst regions. The SED from 
the UV to submm can represent the bolometric luminosity. 



irrespective of the dust extinction. We construct an 'evolu- 
tionary' SED model for starburst regions from a simple but 
realistic point of view. We show that our SED model can ex- 
plain a wide variety of SEDs from optically thin starbursts 
(UVSBGs) to optically thick ones (ULIRGs). It is impor- 
tant to note that our SED model can derive not only the 
SFR, but also the other starburst properties like the optical 
depth. Moreover, we can derive the apparent effective radius 
of starbursts at various wavelengths for a given geometry, 
since we properly take into account the radiative transfer in 
a dusty medium. Therefore, the evolutionary stages of star- 
burst galaxies can be investigated in a unified manner for a 
whole range of activity, which can then be compared with 
the observations of the SED and the efi^ective radius, directly. 
The systematic estimation of the SFRs and the effective ra- 
dius can give an answer to the question: is the panchromatic 
limit of starburst intensity reported by Meurer et al. (1997) 
also valid for ULIRGs? 

This paper is organized as follows. In Section 2, we de- 
scribe our evolutionary SED model for starburst galaxies. In 
Section 3, we summarize the model properties. In Section 4, 
we apply our model to a sample of UVSBGs and ULIRGs 
in the local Universe and describe the resulting properties 
of nearby starburst galaxies derived from our SED model 
fitting. Sections 5 and 6 give discussions and conclusions, 
respectively. 



2 EVOLUTIONARY SED MODEL FOR 
STARBURST REGIONS 

2.1 Star formation and chemical evolution 

We consider a model for the chemical enrichment of a star- 
burst region into which gas is supplied continuously. When 
a large amount of gas is supplied into the central region 
of a galaxy, a starburst is triggered. Since the details of the 
gas supplying mechanism arc yet to be investigated, it would 
not be outrageous to suppose that star formation and chem- 
ical evolution in the starburst regions are approximately de- 
scribed by an infall model of chemical evolution (Arimoto, et 
al. 1992). Ubiquitous galactic outflows from the starbursts 
have been observed with the X-ray emission from the hot 
gas driving the flow, optical line emission produced by the 
warm gas, and the interstellar absorption lines (e.g. Heck- 
man et al. 2000; Lehnert & Heckman 1996). Heckman et 
al. (2000) suggested that most of the outflow gas consists 
of ambient material which has been 'mass-loaded' into the 
hot gas. Therefore, we Eissume that the amount of gas in 
the outflow from a starburst region has negligible impact on 
the chemical evolution as a whole. This is true if the mass 
loss rate from a starburst region is less than 10% of the 
SFR. Under this simple evolutionary picture, a starburst is 
characterized by the rates of gas infall and star formation. 

The time variation of gas mass Mg (t) , total stellar mass 
M«(t), and gas metallicity Zg{t) are given by the equations: 



dt 
dM^{t) 

dt 

d{ZgMg) 
dt 



-i;{t) + E{t) + (1) 

m - m, (2) 

-ij{t)Zg{t) + Ez{t) + ii{t)Zi, (3) 
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where E(t), Ez(t), and Zi are the ejection rates of the gas 
and the metals from dying stars, and the metalUcity of the 
suppUed gas, respectively. We set the initial condition as 
Mg{0) = M«(0) = for all the calculations in this paper. 
Although some amount of gas and stars unrelated to the 
starbursts are initially expected in the starburst region, we 
assume that the amount of initial gas and stars are neg- 
ligible in both the chemical and photometric evolution of 
starburst region. For Zi, the star formation history before 
the starburst event is important. Wc will explicitly note the 
adopted values of Zi later. We assume Zi is constant during 
the starburst event. 

The SFR is given by 



(4) 



where t* and Mg(t) are the star formation time-scaJe and 
the gas mass in the starburst region, respectively. The gas 
supply rate is given by: 



^ / N Mt ft 
5,(,) = _exp 



(5) 



where Mr and U are the initial gas mass in the reservoir 
surrounding the starburst region and the time-scale of gas 
supply, respectively (Arimoto ct al. 1992). Physically, the 
time-scales ti and t, can be expressed in terms of the dynam- 
ical time, sound-crossing time, and cooling time, depending 
on what triggers the starburst. It is however difficult to as- 
certain the characteristic time-scale of each starburst from 
observations. Thus, we hereafter analyse the simplest case, 
in which a starburst is characterized by only one evolution- 
ary time-scale to, i.e., we adopt to = U = tt. 

For all cases in this paper, we adopt the Salpeter ini- 
tial mass function (IMF) with the lower and upper mass 
limit of 0.1 M0 and 60 M0, respectively. Equations (l)-(3) 
are numerically solved by using an evolutionary population 
synthesis code of Kodama & Arimoto (1997). The adopted 
stellar libraries and evolutionary tracks are the same as those 
in the original code. Details of nucleosynthesis prescription 
is given there. 

The effect of dust on the SEDs is significant and pre- 
dominant in starburst regions. In this paper, we adopt a 
simple model in which the dust-to-metal ratio 5o is constant; 
i.e., Moit) = SoZg{t)Mg{t). As described below (Section 
2.3), we use three types of dust model for the Milky Way 
(MW), Large and Small Magellanic Clouds (LMC; SMC). 
The adopted values of So for MW, LMC and SMC are 0.40, 
0.55, 0.75, respectively (Takagi 2001). Starbursts could be 
the most ideal case for constant S{t), since 1) only type II 
supernovae contribute to the supernova rate and 2) the gas 
fraction in molecular clouds is large (see Dwek 1998). 



2.2 Model geometry of starburst region 

We consider a starburst region in which stars and dust are 
distributed within a radius r*. We introduce a mass-radius 
relation for the starburst region; 



n 

Ikpc 



= e 



M, 
IO^Mb 



(6) 



Q. Starburst galaxies arc characterized by a large surface 
brightness, which are roughly constant (Armus, Heckman & 
Miley 1990; Meurer et al. 1995, 1997). In the adopted mass- 
radius relation, 7 = 1/2 gives the constant surface brightness 
for constant O. Therefore, we adopt 7 = 1/2 throughout this 
paper. A variation in the surface brightness can be expressed 
by the different values of in our model. Note that not only 
the surface brightness but also the SED feature is preserved 
for different values of M* if 7 = 1/2, since the source func- 
tion within the starburst is conserved. Therefore, multiple 
systems with the total stellar mass of M,, in which each 
component has the same surface brightness, have the same 
SED as that of a unit system with a stellar mass of M*. 

We assume that the stellar density distribution p(r) is 
given by a generalized King profile; 



p(r) = po 



1 + 



(7) 



where is a compactness factor which expresses the matter 
concentration; the mean density becomes higher for smaller 



where po is the stellar density at the centre of starburst 
region and rc is a core radius of stars. Wo adopt the stel- 
lar density distribution of typical elliptical galaxies; i.e., the 
concentration parameter c = log(rt/rc) = 2.2, and /? = f 
(Combes et al. 1995). Recently, Scoville et al. (2000) and 
James et al. (1999) presented K-hand images of starburst 
galaxies, and demonstrated that luminosity profiles of star- 
bursts are well represented by the r^^* profile which is quite 
similar to those of elliptical galaxies. However, in order to 
derive the stellar density distribution in starbursts, it is im- 
portant to eliminate the efl^ects of dust extinction, since the 
effect of radiative transfer can easily alter the luminosity 
profile from the original one. As we will show later, this 
effect is not negligible even in the /sT-band in ULIRGs (see 
Section 3.3). Therefore, the true density distribution of stars 
in starburst galaxies is difficult to determine directly, and we 
therefore adopt the typical density profile of elliptical galax- 
ies as a first guess. This topic is further discussed in Section 
3.3 and 4.4. 

Following Takagi ct al. (1999), we assume that dust is 
distributed homogeneously within a radius rt of the star- 
burst region. It is plausible that the dust distributes more 
diffusely than the stars, due to feedback from supernovae. 
When the light from centrally concentrated stars dominates 
the SED, i.e., optically thin case, this geometry results in a 
similar SED to the case of shell geometry, which is suggested 
for UVSBGs by Gordon et al. (1997), Meurer, Heckman & 
Calzetti (1999), Witt & Gordon (2000). However, in the op- 
tically thick case, like ULIRGs, the shell geometry results in 
the spectral cut-off around NIR (see the results by Rowan- 
Robinson & Efstathiou 1993), and therefore UV-NIR SED 
of ULIRGs cannot be reproduced with this geometry with- 
out invoking the other components, like underlying stellar 
populations and/or AGNs. As we show in Section 4.2, our 
model can essentially reproduce UV-NIR SEDs of ULIRGs 
only with starburst stellar populations. These SED fitting 
results are confirmed with emission line measurements and 
the observed effective radii in Section 4.3 and 4.4, respec- 
tively. Therefore, we believe that the adopted geometry is 
suitable to approximate the real geometry of starburst re- 
gions with various optical depths. 
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2.3 Dust model 

The dust model is adopted from Takagi (2001) who suc- 
cessfully reproduced the extinction curves observed in MW, 
LMC and SMC, as well as the spectrum of Galactic cirrus. 
The difference among the MW, LMC and SMC extinction 
curves is attributed to the variation of the ratio of carbona- 
ceous dust (graphite and PAHs) to silicate grains. 

The extinction curve is given by the cross-section per 
hydrogen: 

where nn is the number density of hydrogen, crx^k{0') is the 
extinction-cross section of dust particle with the dust con- 
stituent k and size a. The size distribution of each dust con- 
stituent duk/da is taken from Takagi (2001). For a constant 
dust-to-metal ratio, the extinction curve at starburst age t 
is given by: 

(t) = -^o-oA, (9) 

where Z^(t) is the gas metallicity with respect to hydro- 
gen mass (as opposed to the total gas mass) and is the 
metallicity of the ISM under the same definition. We adopt 
Z,'/ = 0.024, 0.011, 0.0034 for the MW, LMC and SMC 
extinction curves, respectively (Pagel 1997; Lequeux et al. 
1979). The time variation of optical depth Tv{t) is given by 
nH{t)(7^{t)rt{t). Note that rt{t) is defined by the time de- 
pendent total stellar mass M,{t). According to the adopted 
mass-radius relation, tv oc nnrt oc Mnr^^ oc /hMtt^^ oc 
fi-iMTifstarMr)-^ oc fHfrtlr wherc Mh and fn are the 
total mass of hydrogen and the mass fraction of hydrogen, 
respectively; therefore, tv does not depend on the value of 
Mt. 



2.4 Intrinsic SED and radiative transfer 

We calculate unobscured stellar SEDs by using the popula- 
tion synthesis code of Kodama & Arimoto (1997), in which 
the effect of stellar metallicity is fully taken into account. We 
solve the equation of radiative transfer by using a computa- 
tional code developed by Takagi (2001). Isotropic multiple 
scattering is assumed and the self-absorption of re-emitted 
energy from dust is fully taken into account. The temper- 
ature fluctuation of very small dust particles is calculated 
consistently with the radiative transfer. 

Wc assume no gradient of the stellar population along 
the radius of starburst region. Although gas emission is not 
considered in our model, a modification of the total SED due 
to gas emission is not significant unless a starburst is con- 
siderably younger than 10 Myr (Leitherer & Heckman 1995; 
Fioc & Rocca-Volmerange 1997). The contribution from gas 
emission to the continuum light is especially important in 
the NIR. As we will show later, no systematic discrepancies 
in the NIR flux are found between model results and obser- 
vations. Therefore, we believe that starburst galaxies are old 
enough to have negligible contribution from gas emission to 
the continuum. 
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Figure 1. Chemical evolution for to =10, 20, 50, 100 and 200 
Myr with Zi=0. a) Evolution of gas mass fraction fgas, stellar 
mass fraction fstar, and metal mass fraction f metal- Both fgas 
and fstar are almost independent of to so that the variation of 
both quantities is no more than the thickness of dashed curves. 
Luminosity-weighted mean stellar metallicities are indicated by 
dot-dashed lines, b) Evolution of optical depth at V-band. For 
two dotted lines, the different values of Zi (O.OIZq for lower and 
O.IZq for upper) are adopted with to = 100 Myr. 



3 MODEL PROPERTIES 
3.1 Evolutionary properties 

Fig. 1 shows how chemical properties and tv evolve as 
a function of t/to, for to = 10,20,50,100, and 200 Myr. 

The gas mass fraction fgas, stellar mass fraction fstar, and 
metal mass fraction f metal are deflned as Mg/Mr, M^/Mt, 
and ZgMg / Mt , respectively. The luminosity- weighted mean 
metallicity is expressed as (Arimoto & Yoshii 1987): 

7 ^_ ^ZstarLstar 

< ^star >— ) (lOj 

^■t-'star 

where Zstar and Lgtar are the metallicity and luminosity 
of each star, respectively. In this section, we simply adopt 
Zi = 0, unless otherwise explicitly noted. 

The variation of fgas and fstar as a function of t/to 
is almost independent of to- After t/to — 2, the evolution 
of fmetai{= fgasZg) and Tv depend little on to. Dust mass 
evolves in proportion to f metal, but ry evolves somewhat 
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differently, since rt increases gradually as stellar mass accu- 
mulates. The effect of Zi on tv is negligible for t/to > 1. 

Fig. 2 shows the evolution of Lboi, Ljr, ijj(t)/Li,oi, 
Lfir/Lh, B — H, r{?, and T"^ . is the effective optical 
depth at V-band defined as Lv = i'v exp(— r{J?), where 
and Lv are the intrinsic and attenuated luminosity, respec- 
tively, is the effective dust temperature which gives the 
same peak wavelength of dust emission as that of a model 
SED with an emissivity law oc A~^. 

We summarize the model properties as follows: 

1) The IR luminosity evolution suggests that starburst 
galaxies should be found mostly in the late phase (i/io ^ 2), 
after the SFR and the dust mass become maximum. The 
asymptotic relation between SFR ^{t) and L^oi for the late 
phase t/to > 2 is 

Lboi ^ ^(t) t 
1.7 X IO^Lq IMoyr-i to' ^ ' 

We also derive an asymptotic relation for the mass-to-light 
ratio for the late phase t/to > 2; 

i^=8.05xlO-(f)^Y-^V"Mo. (12) 
Lboi \toJ l^lOOMyry Lq ^ ' 

These ratios will be used in the analysis of starbursts later 
(Section 5). 

2) The variation of Lfir/Lh as a function of t/to is 
nearly independent oito- Lfir/Lh decreases monotonically 
with time. Due to this weak dependence of SED quantities 
on to, it is difficult to determine the absolute starburst age 
from the SED fitting alone. 

3) The dependence of B — H on to is similar to that of 
TV and . 

Fig. 3 shows the SED evolution for to = 100 Myr, 9 = 1, 
the MW extinction curve, and Zi = 0. The SEDs arc quite 
similar to each other at t/to < 2. At t/to > 2, the peak 
wavelength of dust emission increases and the FIR excess 
(Lfir/Lh) decreases. 



3.2 Dependence on compactness factor 

Fig. 4 shows how the characteristic SED quantities depend 
on the compactness factor 0. While Lfir is not sensitive to 
0, Lfir/Lh increases significantly for < 1. At a small 
(i.e., optically thick), the SED at optical wavelengths asymp- 
totically reaches the source function of the stellar light, since 
the light mainly comes from near the surface for such a 
highly obscured region. This is the reason why a starburst 
region becomes bluer for O<0.5. 

There are two possible reasons for large values of tv; 

1) starbursts occur in metal rich clouds (Zi > 0.1.^©) and 
are still very young {t/to < 0.3) as shown in Fig. lb, and 

2) geometry of starbursts is compact (0 < 1) as shown in 
Fig. 4d. If the former reason is the case, ULlRGs should 
be younger than UVSBGs, systematically. However, we find 
no such systematic difference as we show in Section 4. In- 
stead, SED fitting results in smaller for ULIRGs than 
those of UVSBGs. Thus, the variation of tv in starbursts is 
mainly caused by the geometrical effect; i.e., O. This result 
is reasonable, since the luminosity evolution suggests that 
starbursts are less luminous in the early phase of starbursts. 
In later phases, the high luminosity and large fraction of 



life time during this phase as a starburst would increase the 
chance of the starbursts being observed. The importance of 
on the variation in starburst SED is illustrated clearly in 
Fig. 5. 

3.3 Effective radius in various wavelengths 

Fig. 6 shows the growth curve of the apparent surface bright- 
ness at 0.22, 0.7, 1.1, 2.2, and 100 /xm for the models of UVS- 
BGs {t/to = 2.0, = 1.5 and the MW extinction curve) 
and ULIRGs {t/to = 2.0, = 0.4 and the SMC extinction 
curve). Here we adopt the typical values of model parame- 
ters for each UVSBGs and ULIRGs, which arc derived from 
the SED fitting described in the next section. At most of 
the observed wavelengths, as seen in Fig. 6, the apparent 
effective radii arc systematically larger than the intrinsic ef- 
fective radius which corresponds to the case without any 
extinction. In the FUV to NIR wavelength region in which 
the dust emission is negligible, the apparent effective radius 
increases with decreasing wavelength. This is because the 
contribution from stars in the central region decreases with 
increasing optical depth toward the shorter wavelength re- 
gion. Therefore, only the growth curve in 2.2 yum follows that 
of the input stellar distribution for the case of UVSBGs. 
Meurer et al. (1995) presented the surface brightnesses of 
UVSBGs at FUV (0.22 /im) obtained with the HST/Faint 
Object Camera (FOG). These surface brightness is not well 
represented by the r^^'* or exponential profile. This result is 
consistent with the model for UVSBGs in which the growth 
curve at 0.22 pm is significantly different from the intrin- 
sic growth curve (dashed line). Note that both the r^^'* and 
exponential profile have been used for normal gala^xies in 
which the effect of dust is not significant, compared to star- 
burst galaxies. On the other hand, in the case of ULIRGs, 
the apparent effective radius even at 2.2 jim becomes twice 
the intrinsic effective radius, although the growth curve is 
somewhat similar to the intrinsic one. In both UVSBGs and 
ULIRGs, the effect of scattering on the growth curve is neg- 
ligible. The apparent effective radius at 100 /im is also differ- 
ent from the intrinsic effective radius, since it is determined 
by the distribution of dust, rather than that of stars. 



4 SED FITTING TO UVSBGs AND ULIRGs 

We apply our evolutionary SED model to two samples of 
local starburst galaxies; UVSBGs and ULIRGs. All models 
are calculated with to = 100 Myr, Z, =0.1 Zq. A definitive 
value of to is only important for the absolute time scale of 
starburst events which is not discussed here. The value of 
Zi should vary for each starburst galaxy, since Zi depends 
on the star formation history in precursors of a starburst 
galaxy. We expect that the gas fraction in precursors of star- 
bursts should be large, and therefore assume that the level 
of chemical enrichment in precursors is low and comparable 
to that in the SMC. Our conclusions are independent of the 
choice of Zi. 

Fitting parameters are the starburst age t, and the 
compactness factor 0. We choose the best-fitting extinction 
curve from the MW, LMC and SMC type. The SEDs are 
normalized by the initial mass Mr of a gas reservoir. The 
best-fitting model for each galaxy is sought by a standard 
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Figure 2. The evolution of characteristic quantities of starburst galaxies, a) -Lfcoi (solid line) and Liji (dashed line), b) the ratio of SFR 
to Ljoj, c) Lpi f! / Lh ^ d) B — H, e) , f) effective temperature of dust Tf^ for various to (open circles for to= 10 Myr, squares for 20 
Myr, triangles for 50 Myr, solid circles for 100 Myr, and diamonds for 200 Myr). L/^ and Lpjn are the luminosity in the wavelength 
range of 8 - 1000 fira and 42 - 122 fj,ra, respectively. We adopt the MW extinction curve, = 1, Mt = 1 Mq and Z, = 0.0. 



least square method. We adopt the cosmology of Ho = 75 
km sec~^ Mpc~^ and qo = 0.5 throughout this paper. 



4.1 UVSBGs 

Wc analyse the same sample of UVSBGs analysed by Tak- 
agi et al. (1999). Photometric data from FUV to NIR are 
taken from Gordon ct al. (1997) who presented a homo- 
geneous sample of aperture-matched SEDs of 30 starburst 
galaxies. All of the sample was detected by the /[/i? satellite. 
UV-sclcctcd starbursts arc relatively bright in the UV and 
expected to be less attenuated by dust. The aperture size 
of the FUV-NIR photometry is 10" x 20", which is equal to 
~ 4.5 kpc at their median distance of 60 Mpc. The IRAS flux 



densities are taken from NED (NASA/IPAC Extragalactic 

Database). The IRAS fluxes are derived for the total galaxy, 
while the UV-NIR photometric data arc aperture limited in 
the central region of starburst galaxies on average. Never- 
theless, the aperture correction does not change our main 
conclusions, since the FIR emissions tend to concentrate at 
galactic centre in the galaxies analysed here (Calzetti et al. 
1995). We fit the SEDs of 16 galaxies out of 30, which were 
detected in all four IRAS bands (12/Ltm, 25/^m, 60/Ltm and 

First, the best-fitting model is determined by the SED 
fitting from the U-hand to the FIR, i.e., without FUV data 
obtained with the lUE satellite. The fitting results from the 
?7-band to the FIR are shown in Fig. 7 with dashed lines. 
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Figure 3. The spectral evolution of a starburst for = 1 with Mt=1 M©, the MW extinction curve and Zi = 0. The evolutionary 
stage for each line is denoted near the lines. For clarity, the SEDs are shifted vertically with the denoted value in parenthesis at the 
emission peak. 



All the best-fitting models result either in reasonable fits or 
in underestimates to the FUV data. 

This discrepancy in the FUV flux can be reasonably 
explained by the effect of photon leakage. Although we con- 
sider a homogeneous distribution of dust, regions with low 
gas density are naturally expected in the starbursts, because 
of a cavity produced by successive supernova explosions, for 
example. The effect of photon leakage is expected to be more 
significant in the FUV compared with the longer wavelength 
region, since the starburst population is intrinsically lumi- 
nous in the FUV. Therefore, the fitting results are consistent 
with our proposed picture of the starburst region. 

We estimate the effect of photon leakage as follows. 
Although various levels of extinction for leaking photons 
are expected, depending on the degree of inhomogeneity in 
the dust distribution, the contribution of leaking photons 
from optically thiner regions could be much more signifi- 
cant. Therefore, for simplicity we assume that the extinc- 
tion of leaking photons is negligible. We define the energy 
fraction of photons which leak out from the stars directly 
through the starburst region, /loak, as: 



n = (1 - /leak)i^A(T-A) + fle^kF^irx = 0), 



(13) 



where F^*"* is the observed fiux, and Fx{t\) is the model 
flux for a given optical depth without considering the photon 
leakage. By using F\{tx) and Fx{tx = 0) for each starburst 
galaxy derived from the (7-band to FIR SED (i.e., without 
FUV data), we determine the most plausible value of /leak 
by the least square method by using all photometric data 
including the FUV; i.e., in this fitting process, wo use /loak 
as a free parameter. Also, Mt is re-adjusted to give the least 
square value for the SED with the leakage correction, which 
is necessary to derive a more accurate bolometric luminosity, 
and therefore SFR. 

The best-fitting models with the photon leakage are in- 
dicated by solid lines in Fig. 7. As a whole, these models 
show a good agreement with the observed SEDs. For most 
of the sample galaxies, the SEDs in the FUV can bo repro- 
duced with a /leak less than 0.1. Note that the correction 
of SED due to the photon leakage effects only the FUV. In 
Table 1, the best-fitting parameters are given together with 
the resulting Ltoi, M*, Mo, SFR and rf. From the SED 
fittings with our adopted evolutionary time-scale to = 100 
Myr, most of the derived starburst ages are within 200 — 400 
Myr. As we already noticed, starbursts are most luminous 
during the ages of t/to ~ 2 — 4. Therefore, the derived ages 
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Figure 4. The dependence of a) Lpm (thin lines) and Lm (thick lines), b) Lfir/ Lh , c) B — H, d), -ry (thin lines) and (thick 
Unes), e) on 0. The case for to = 100 Myr are shown in each panel. Sold, dashed and dot-dashed lines indicate the case of t/to = 1-0, 
3.0 and 5.0, respectively. We adopt the MW extinction curve, Mt = 1 M©, and Zj = 0.0. 



with the SED fittings are reasonable for the staxbursts. For 
all the sample but three, we found that the MW extinction 
curve is most suitable for UVSBGs. Note that the absorp- 
tion feature at 2175A is not used to determine the extinc- 
tion curve. The MIR colours are important for the resulting 
extinction curve, since the difference in the relative frac- 
tion of dust constituents can effect the MIR SED in which 
very small graphite grains and PAHs play significant role. 
The mean value of Q is 1.6, corresponding to = 1.2 for 
t/to = 2.0 in the case of the MW extinction curve. 

Some galaxies, such as NGC5236, NGC6052 and 
NGC6217, show a notable excess of 100 /im flux compared 
to the model SEDs. As we noted, the aperture size adopted 
for IRAS data is larger than those of the FUV - NIR data. 
Therefore, it is likely that the FIR emission from cold dust. 



heated by the interstellar radiation field generated by the 
underlying stellar populations, contributes to the 100 /xm 
flux. 

For IC214, the model with the photon leakage still un- 
derestimates the FUV fluxes. Since this galaxy along with 
NGC1614 is the faintest in FUV among the UVSBG sample, 
the signal-to-noise ratio of the lUE observation is low (3.2 
for 1515A; Kinney et al. 1993). A more sensitive observation 
would be necessary for a detailed discussion on the origin of 
this discrepancy. 

4.2 ULIRGs 

The sample of galaxies with Lir > 10^^ Lq, i.e., ULIRGs, 
analysed here were originally studied by Sanders et al. 
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Figure 6. Growth curves of starburst region at various wavelengths. The vertical axis is the fraction of flux emitted within the projected 
radius r. Dashed line is for the case without dust. Effective radii for each wavelength is indicated by solid circles. Dotted line indicates 
the growth curve at 0.22 /im without a contribution from scattering light. Left panel: the case for UVSBGs (t/to = 2.0, = 1.5 and the 
MW extinction curve), right panel: the case for ULIRGs (t/to = 2.0, = 0.4 and the SMC extinction curve). 
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Figure 7. Fitting results for UV-sclcctcd starburst galaxies. The dashed lines indicate the fitting results without taking into account 
the photon leakage. In this case, we fit the model SEDs to observations, except for FUV data. The solid lines cire the best-fitting models 
for observed SEDs including the FUV data, considering the photon leakage. XLx is in solar unit. 



(1988). Although a starburst is the only energy source in 
our model, wc also apply it to Scyfcrt galaxies in which the 
AGN could considerably contribute to the dust heating. Ac- 
cording to Veilleux et al. (1999a), IRAS05189, Mrk231 and 
Mrk273 are Seyfert galaxies, although the line-to-continuum 
ratio of the 7.7 /xm emission feature indicates that Mrk273 
is a starburst-dominanted ULIRG (Rigopoulou et al. 1999). 

All photometric data from the B-band to the submm 
wavelengths are taken from Rigopoulou, Lawrence & 
Rowan- Robinson (1996). The adopted aperture size is 5" for 
all ULIRGs, corresponding to ^ 5 kpc for the mean redshift. 
The adopted aperture is the smallest among those tabulated 
in Sanders et al. (1988), where the original photometric data 



are taken from. Although IRAS and submm data are ob- 
tained with larger apertures, the compactness of emission 
source ensures negligible corrections for the aperture mis- 
match (Sanders et al. 1988; Wynn- Williams & Beckhn 1993; 
Soifor et al. 2000). 

For ULIRGs, galaxies with secure FUV photometry are 
relatively rare. Thus, we simply determine the best-fitting 
model for all the photometric data available. We find that 
our models again give a flux lower than that observed in the 
shorter wavelength region. We estimate /loak for ULIRGs by 
attributing this deficiency to the photon leakage. 

The best-fitting models without considering the leak- 
age are shown in Fig. 8 with dashed lines. Our model can 
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Figure 7. - continued. 



reproduce the observed SEDs especially well at wavelengths 
longer than NIR, except for IRAS05189 and Mrk231, which 
are candidates for AGN powered ULIRGs. For both galaxies, 
it is found that the observed flux from NIR to MIR show 
clear excess comparing with the best-fitting model, which 
could be accounted for by the emission from hot dust around 
AGN. Thus, our fitting results arc consistent with the spec- 
troscopic classification of ULIRGs. All of starburst powered 
ULIRGs arc best-fit with the SMC type extinction curve 
except for IRAS08572. This shows a clear contrast to the 
results for UVSBGs, in which the MW extinction curve is 
found to be better. All the best-fitting models arc character- 
ized by the small compactness factor O ~ 0.5, corresponding 
to TV ~ 20; as a result, the strong silicate absorption feature 



at 9.7 /um is a common factor in ULIRGs. This result sug- 
gests that the starburst regions of ULIRGs are much more 
compact in size compared with the UVSBGs. 

As already mentioned above, we can see that the best- 
fitting models tend to give a flux in the optical region lower 
than the observed one. Following our treatment of UVSBGs, 
we determine /icak with the least square method, assuming 
negligible extinction for leaking photons. The best-fitting 
models with the leakage are shown in Fig. 8 with solid lines. 
The resulting fitting parameters are given in Table 2. Due 
to the simple assumption that the extinction is negligible 
for the leaking photon, the resulting SEDs in the UV are 
significantly bluer. Although the values of /leak are smaller 
than those for UVSBGs, the effect of the leakage on the 
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Figure 7. - continued. 



Table 1. Fitting results for UVSBGs 



name 


t/to 


e 


EC 


/leak 


Mt 


Lboi 


M» 


Md 


SFR 














[Mq] 


[Lq] 


[Mo] 




[M0 yr-i] 




IC214 


4.0 


0.8 


MW 


0.012 


2.2 10^" 


2.1 10" 


2.0 lO^" 


3.8 10^ 


25.2 


2.2 


NGC1140 


3.0 


1.4 


LMC 


0.165 


3.5 10* 


4.9 10'^ 


2.7 10* 


1.1 10^ 


0.7 


0.4 


NGC1569 


0.7 


2.6 


SMC 


0.101 


8.9 lO'^ 


1.4 10^ 


1.3 lO'^ 


1.7 10^ 


0.3 


0.3 


NGC1614 


2.0 


0.8 


LMC 


0.002 


1.8 IQi" 


3.5 10" 


1.0 lO^" 


6.1 lO'' 


57.7 


2.5 


NGC4194 


2.0 


1.5 


MW 


0.004 


3.7 10^ 


7.1 IQio 


2.1 10^ 


9.1 10^ 


11.7 


1.6 


NGC4385 


3.0 


1.3 


MW 


0.054 


1.1 10^ 


1.5 IQio 


8.1 10* 


2.4 10® 


2.1 


1.3 


NGC5236 


2.0 


1.6 


MW 


0.028 


1.2 lO'' 


2.4 10i° 


6.9 10* 


3.0 10*^ 


3.9 


1.2 


NGC5253 


0.3 


2.5 


MW 


0.084 


2.3 10** 


1.4 10^ 


7.9 10^ 


7.3 10-* 


0.5 


0.9 


NGC5860 


4.0 


1.3 


MW 


0.035 


3.6 10'^ 


3.4 10^" 


3.2 10" 


6.2 10® 


4.1 


0.9 


NGC6052 


2.0 


1.8 


MW 


0.038 


3.2 10'^ 


6.0 10^" 


1.8 10" 


7.7 10® 


10.0 


0.9 


NGC6090 


2.0 


1.5 


MW 


0.017 


1.4 10^" 


2.7 10" 


7.9 10" 


3.5 10^ 


45.0 


1.5 


NGC6217 


3.0 


1.3 


MW 


0.013 


9.6 10* 


1.3 10^0 


73 10* 


2.2 10® 


1.9 


1.4 


NGC7250 


2.0 


2.4 


MW 


0.088 


1.5 10* 


2.8 10^ 


8.1 10^ 


3.5 10^ 


0.4 


0.5 


NGC7552 


2.0 


1.3 


MW 


< 0.001 


4.6 10^ 


8.8 10^° 


2.6 10^ 


1.1 lO'' 


14.5 


2.1 


NGC7673 


2.0 


2.2 


MW 


0.032 


1.3 10^ 


2.5 IQio 


7.3 10* 


3.2 10® 


4.1 


0.6 


NGC7714 


2.0 


1.8 


MW 


0.057 


2.5 10^ 


4.7 IQi" 


1.4 10^ 


5.9 10® 


7.7 


0.9 



SED is prominent because of heavy obscuration of the stel- 
lar continuum. In order to verify the fitting results, we plot 
the FUV photometric data (Trentham et al. 1999; Goldader 
et al. 2002) for some of the ULIRG sample with open cir- 
cles in Fig. 8. For Mrk273, IRAS08572 and IRAS12112, the 
FUV data taken by the HST arc consistent with the best- 
fitting model without considering the leakage (dashed line) , 



while for IRAS15250 and IRAS22491 the FUV data can be 
explained by the best-fitting model with the leakage (solid 
line). Only for Arp220, the FUV data cannot be explained 
by both cases. According to Goldader et al. (2002), these 
data are problematic due to uncertainty of sky level, since 
the angular size of Arp220 is large enough to cover entire 
field of view of HST/STIS. Thus, fitting results are con- 
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Figure 8. The same as Fig. 7, but for ultraluminous infrared galaxies. Open circles for IRAS12112, IRAS08572, IRAS15250, Arp220 and 
IRAS22491 indicate the photometric data, obtained by the Hubble Space Telescope (Trentham, Kormendy & Sanders 1999; Goldader et 
al. 2002). 



sistent with FUV observations, except for the problematic 
case of Arp220, when we consider the derived value of /leak 
as the upper limit for some galaxies. 

As mentioned earlier, starburst ages are difficult to as- 
sign, but our resulting ages are in the range of 50 — 500 Myr, 
which shows no systematic deviation from those derived for 
UVSBGs. For starburst dominated ULIRGs, the mean value 
of and is 0.46 and 4.9, respectively. Thus, ULIRGs are 
more compact than UVSBGs by a factor of 3.5. 

For most of the ULIRGs, the SMC extinction curve is 
found to provide the best fit. Note that neither the obser- 
vations of silicate absorption features nor the bump at 2175 
A are used to distinguish the type of extinction curve in 



the SED fitting. The SED from MIR to submm wavelengths 
can be better reproduced with the model in which the self- 
absorption of MIR emission by dust is significant; i.e., the 
model with the SMC extinction curve. 



4.3 Comparison with emission line measurements 

4-3.1 Star formation rate 

In Fig. 9, we compare the resulting SFRs with the ob- 
served values derived from Ha luminosities for UVSBGs 
(Storchi-Bergmann et al. 1994) and from the FIR luminos- 
ity of ULIRGs by using the relation presented by Kennicutt 
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Figure 8. - continued. 
Table 2. Fitting results for ULIRGs 



name 


t/to 


e 


EC 


/leak 


Mt 




M* 


Md 


SFR 














[Mq] 


[L0] 


[Mq] 


[M©] 


[Me/yr] 


UGC05101 


5.0 


0.3 


SMC 


0.005 


1.1 10^1 


2.1 10^^ 


1.0 10" 


2.5 10* 


69. 


4.7 


IRAS12112 


2.0 


0.3 


SMC 


0.006 


7.8 


1.5 10^2 


4.3 10"^ 


3.5 10* 


244. 


5.4 


Mrk273 


2.0 


0.4 


SMC 


0.008 


6.6 10^° 


1.3 10^2 


3.7 10^° 


3.0 10* 


207. 


4.9 


IRAS08572 


0.5 


0.8 


MW 


0.036 


9.2 10i° 


1.7 10i2 


7.7 109 


6.0 10^ 


285. 


4.9 


IRAS14348 


2.0 


0.3 


SMC 


0.008 


1.0 10" 


2.0 10^2 


5.8 10^° 


4.8 10* 


327. 


5.4 


IRAS15250 


1.0 


0.6 


SMC 


0.005 


4.8 10^° 


9.1 10" 


1.2 10^° 


1.4 10* 


185. 


4.5 


Arp220 


0.7 


0.4 


SMC 


0.004 


7.2 10i° 


1.4 10^2 


1.0 10^° 


1.4 10* 


260. 


5.4 


IRAS22491 


1.0 


0.6 


SMC 


0.011 


5.7 IQio 


1.1 IQi^ 


1.4 IQlO 


1.6 10* 


221. 


4.4 


IRAS05189" 


4.0 


0.4 


MW 


< 0.001 


1.4 10" 


2.7 10^2 


1.2 10" 


2.4 10* 


158. 


4.7 


Mrk231'" 


2.0 


1.0 


MW 


< 0.001 


1.6 10" 


3.1 10^2 


8.9 10^" 


3.9 10* 


506. 


2.8 



Note: " candidates for AGN powered ULIRGs 



(1998). The SFRs in our model are basically determined 
from the bolometric luminosity from UV to FIR. Therefore, 
the observed SFRs for ULIRGs should be almost identical to 
those from the SED fitting since the emission from the FIR 
dominates the SED of ULIRGs. Thus, we can see clearly this 
agreement at the high SFRs. On the other hand, our SFR 
derivation is independent of that from Ha line emission for 
UVSBGs. Although the difference between the model results 
and observation becomes large for galaxies with the SFR < 
1 Mq yr"\ the resuhing SFRs for UVSBGs agree weU with 
those derived from Ha and IR luminosities over three orders 
of magnitude without any systematic difference. Therefore, 



the SFR estimated by the SED fitting can be applied well 
to any starburst galaxy regardless of their optical depth. 

Our SED model suggests that leaking photons con- 
tribute to the UV luminosities of starbursts. Therefore, if 
the UV luminosity is corrected for dust extinction without 
considering this contribution, the UV-derived SFR can be 
easily overestimated. This effect may explain the excess of 
UV luminosities for a given Hq luminosity, reported by Sul- 
livan et al. (2000). We emphasize that, in order to derive 
SFRs from the continuum emission, bolometric luminosity 
should be used to avoid systematic effects. 
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Figure 9. Comparison of the SFR in MQjrr"^. The observed 
SFRs are derived from the extinction corrected Ha and FIR lumi- 
nosity for UVSBGs and ULIRGs, respectively. AGN candidates, 
IRAS05189 and Mrk231, are not plotted. Solid line indicates the 

linear relation. 
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Figure 10. Comparison of the extinction. The observed colour 
excesses E(B — y)gas are derived from the ratios of emission 
lines. For UVSBGs, Ha/H/3 is used, taken from Storchi-Bcrgmann 
ct al. (1994). For ULIRGs, only galaxies which have observed 
value of Paa are plotted; depicted numbers indicate galaxies; 1 
- IRAS08572, 2 - IRAS12112, 3 - Mrk273, 4 - IRAS15250, 5 - 
IRAS22491. Filled and open squares indicate E{B—V)gae derived 
from Paa /Ha and Ha/H/?, respectively (Veilleux et al. 1999b). 
Solid and dashed line indicate the relation for the ratio of total- 
to-selective extinction Av/E{B-V) = 3.08 for MW and 2.93 for 
SMC (Pei 1992), respectively. The derived relation between 
and E{B — V')gas should be similar to this relation in the case 
that the distribution of ionized gas is the same as that of stars. 



4-3.2 Optical depth and extinction 

The effective optical depth can be determined independently 
of the SFR. In Fig. 10, wo compare the effective optical depth 
at F-band with the colour excess E{B — V)ga.s derived 
from the Balmer line ratios. We see that Ty? is consistent 
with E{B — V)gas- Calzetti et al. (1994) suggested that the 
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Figure 11. Comparison of X-band surface brightness profile of 
NGC1614. Data plotted are taJjen from Alonso-Herrero et al. 
(2001) obtained with the F5T/NICMOS. Solid line is the sur- 
face brightness of the best- fitting SED model for NGC1614. 

extinction inferred from the line ratios in starburst galaxies 
is typically twice as high as that derived from the stellar 
continuum. For optically thick UVSBGs in our sample, the 
resulting tends to be slightly less (within a factor of 2) 
than the expected value for the MW extinction curve and 
therefore, consistent with the results of Calzetti et al. (1994). 

The observed colour excesses E{B — V)gas of ULIRGs 
should be derived from the less attenuated lines, such as 
Paschen and Brackett lines, since the Balmer line ratios 
could be saturated due to the large optical depth. Veilleux 
et al. (1999b) presented E{B — F)gas measurements derived 
from both Paa/Ha and Ha/H/3. Both are plotted in Fig. 
10, connected with dotted lines for each galaxy. Half of our 
sample of galaxies have the observed Paa/Ha. Except for 
IRAS22491, rf of the ULIRG sample are consistent with 
the simply extrapolated relation from the less attenuated 
UVSBGs. For IRAS22491 (the only HII galaxy in our sam- 
ple), E{B — y)gas derived from Paa/Ha and Ha/H,5 differs 
in the opposite sense. Galaxies showing such anomalous be- 
haviour are predominantly Seyfert 2 galaxies in the sample 
studied by Veilleux et al (1999b), and HII galaxies consti- 
tute only a small fraction. The extinction of ULIRGs could 
be investigated more accurately by spectroscopic studies at 
longer wavelengths with higher sensitivity, by future FIR 
satellite missions, such as SIRTF and ASTRO-F. 

In summary, the SED fitting model can reproduce well 
the relation between E{B—V)ga.B and t^?, and therefore, the 
dimming and reddening effects of the light from starburst 
regions are consistently introduced in the model. 

4.4 Comparison of apparent efFective radius 

Since we applied the model not to a whole galaxy, but to 
the central starburst region, the effective radius derived from 
the SED fitting should be compared with that of the central 
starburst region. For the UVSBG sample, unfortunately no 
homogeneous presentation of effective radii of starburst re- 
gions is found in the literature. Prugniel & Heraudeau (1998) 
presented a catalogue of aperture photometry of 7744 galax- 



16 T. Takagi, N. Arimoto and H. Hanami 



6 




f . I . . . . I . . .^~r~~r~— 

1.0 1.5 2.0 2.5 



shown in the Fig. 12, the effective radius of isolated star- 
burst region in the if-band is < 0.7 kpc, which corresponds 
to an intrinsic effective radius (i.e., for the case without dust) 
of ~0.3 kpc. Scoville et al. (2000) have shown that the hght 
profiles of these ULIRGs can be also fitted well with the r^^* 
law, which can be reproduced by the King profile. Thus, the 
stellar distribution in a starburst region cam be roughly rep- 
resented with our model geometry of the King profile with 
log(rt/rc) = 2.2 for ULIRGs. 

As shown by Scoville et al. (2000), many ULIRGs have 
multiple starburst regions - such as Arp 220. As long as each 
starburst region has a constant surface brightness, the SED 
feature is not affected by this multiplicity as remarked in the 
model description (Section 2.2). Thus, the model geometry 
can still work as long as the stellar distribution of individual 
starburst regions are considered. 



Figure 12. Comparison of the apparent effective radii at 1.1/um, 
1.6/im and 2.2/tm with the observations (Scoville et al. 2000; 
solid circles) for galaxies in which the effective radius is deter- 
mined for a single nucleus (see text in detail). For clarity, we add 
the value noted in parenthesis to the effective radii. 



ies and derived _B-band effective radii for a whole galaxy, in- 
cluding all our sample UVSBGs, but three. Note that these 
effective radii cannot be compared with the model results di- 
rectly, since the underlying population could enlarge the B- 
band effective radii, systematically. It is found that B-band 
effective radii published by Prugniel & Heraudeau (1998) are 
systematically larger than those derived from the SED fit- 
ting by a factor of three. However, it is not yet clear whether 
the contribution from underlying population is enough to 
explain this discrepancy or not. This discrepancy may be 
partly due to the adopted concentration parameter c being 
too large for UVSBGs. Alonso-Herrero et al. (2001) studied 
NGC1614 in detail with HST/NICMOS observations, and 
presented a surface brightness profile of the central star- 
burst region in the i^-band. We compare the Jf-band surface 
brightness profile of NGC1614 with that of the best-fitting 
SED model in Fig. 11. Clearly, the model result is in good 
agreement with the observation. Since this galaxy is rela- 
tively luminous and obscured for a typical UVSBG, wo need 
a more comprehensive study using images with high spatial 
resolution to confirm the stellar distribution of UVSBGs un- 
ambiguously. 

Scoville et al. (2000) presented NIR images of ULIRGs 
obtained with the HST/NICMOS, and derived the effective 
radii at 1.1/im, 1.6/im and 2.2/xm. In most cases, the ob- 
served effective radii are less than 1", and are therefore con- 
sistent with the adopted aperture size as starburst regions 
(5") in the SED fitting. We focus on the ULIRGs for which 
the effective radii are determined for isolated major star- 
burst regions, since the multiple structure could enlarge the 
effective radius, systematically. In Fig. 12, the effective radii 
of the models for ULIRGs are compared with the apparent 
effective radii from the HST/NICMOS images. As a result, 
we can see that the radii from the models and the observa- 
tions are consistent with each other within a factor of 2. As 



4.5 Correlations between starburst properties: 
signature of bimodality 

In Fig. 13a, the SFR is plotted against and tv ■ Clearly, 
we can see a power-law correlation between SFR and rJJ?; 
SFR oc {rf f. A similar trend has already boon shown by 
Heckman et al. (1998) and Adelberger & Steidel (2000) be- 
tween Lfir/I/uv and Luv+fir, which are related to the red- 
dening effects and SFR respectively, although the dispersion 
of the correlation they presented was much larger than that 
shown in Fig. 13a. This trend seems to suggest the Schmidt 
law for star formation, in which SFR increases with increas- 
ing gas density. However, note that not t^, but tv is pro- 
portional to the gas column density in the starburst regions. 
The discrepancy between and tv can be clearly seen at 
high optical depth > 5; as a result, the SFR seems to sat- 
urate at TV ^ 5. This suggests that the mode of starburst 
activity changes at a critical column density corresponding 
to TV — 5. 

Such bimodal feature can be seen in the relation be- 
tween the SFR and the compactness factor as shown in Fig. 
13b, in which we can also find a critical compactness fac- 
tor ~ 3. Both bimodal features of the SFR against the 
column density and the compactness suggest that the gas 
density in the starburst regions can directly affect the star- 
burst activity since in general the geometrical compactness 
with the high column density implies the existence of a high 
density region. 

In Fig. 13c, wo plot /icak against Ty? and tv- There 
is an anti-correlation between these properties for tv < 5, 
while this anti-correlation breaks for larger tv. Feedback 
in starbursts with low ISM density, such as UVSBGs, may 
naturally generate regions of porous ISM that provide path- 
ways for UV photon leakage even from central region of star- 
bursts. On the other hand, for starbursts with a high density 
ISM such as ULIRGs, leakage can be expected only from the 
outermost surface of a starburst region and may have little 
to do with the main starburst (Goldader et al. 2002). Thus, 
the break of anti-correlation between /leak and tv can be 
explained by the feedback effects on different ISM density. 
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Figure 13. Relations between the derived quantities. Circles and squares indicate UVSBGs and ULIRGs, respectively, a) SFRs are 
plotted against the effective V-band optical depth (open symbols) and y-band optical depth (solid symbols). The solid line is a power- 
law relation, SFR oc (t^)^. b) SFRs are plotted against 0~^, which is proportional to Ty for constant dust mass, c) Leaking fraction 
as a function of Ty? (open symbols) and tv (solid symbols) . Upper limit on /leak is indicated with an arrow. 



5 DISCUSSION 

In this section, wo will discuss the bimodality found in the 
correlations between starburst properties with our SED di- 
agnostics. First, we focus on the correlation between the 
intrinsic surface brightness and the size of starburst region, 
which shows the bimodality with a clear limit on the size 
of starburst regions. Then, wc investigate the origin of the 
bimodality in the starburst activity by considering the phys- 
ical processes of starburst regions, such as gravitational in- 
stability of the disc and feedback effects. 



5.1 Bimodal starburst intensity 

In Fig. 14, the bolometric luminosities of starburst galax- 
ies are plotted against the effective radii. Filled circles and 
squares are bolometric luminosities and intrinsic effective 
radii (i.e., for the case without dust) derived from the SED 
model of UVSBGs and ULIRGs, respectively. The bolomet- 
ric surface brightnesses of UVSBGs arc found to bo almost 
constant at lO^'^Lokpc"^. On the other hand, a clear limit 
on the intrinsic effective radius can be seen at ~ 0.3 kpc. 
Moreover, the bolometric surface brightnesses of ULIRGs 



are an order of magnitude higher than those of UVSBGs. 
This suggests that the mode of starbursts could change at 
Lboi — 10^^ Lq and Ve ^ 0.3 kpc. 

A limiting bolometric starburst intensity L^oir^^ = 
2 X 10" Lokpc"^ was proposed by Meurer et al. (1997). 
From this intensity limit, they also proposed a global star 
formation limit ipr^'^ ~ 45M0kpc~^yr~^. In Fig. 14, all 
open symbols other than open squares indicate the origi- 
nal data used by Meurer et al. (1997). For open squares, 
we adopt the effective radii of ULIRGs observed with the 
FSr/NICMOS at 2.2/im taken from Scoville et al. (2000). 
These effective radii are systematically smaller than those 
derived from il-band images (Lehnert & Heckman 1996; 
small open diamonds) and from images of Ha emission (Ar- 
mus et al. 1990; small open triangles), although the similar 
class of galaxies is investigated. Thus, from the observational 
results it is suggested that the observed surface brightness 
depends on the observed wavelength. For UVSBGs, the bolo- 
metric surface brightness calculated with the intrinsic effec- 
tive radius is higher by a factor of 5 than the proposed value 
by Meurer et al. (1997). The difference between the intrin- 
sic and observed surface brightness is more remarkable for 
ULIRGs and indicates a bimodal starburst activity. 
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Figure 14. The bolomctric luminosities of starburst galaxies plotted against effective radii. Filled circles and squares are the bolometric 
luminosities and intrinsic effective radii derived from the SED of UVSBGs and ULIRGs, respectively. Open symbols indicate the observed 
bolometric luminosities and apparent efEective radii at 0.22/im (open circles; Meurer et al. 1995), 2.2/im (squares; Scoville et al. 2000). 
Small diamonds indicate the ii-band effective radii and Lh^i taken from Lehnert & Heckman (1996) in which Lj,o; is estimated from FIR 
luminosity (Lf,o( = l.bLpjji; Meurer et al. 1997). When i?-band data is not available, we convert the effective radius measured with 
Ha imaging into R-band effective radius with the average ratio of Ha to R-band effective radius, 0.6 according to Lehnert & Heckman 
(1996). By using the same ratio of effective radii, Ha effective radii taken from Armus et al. (1990) are converted to _R-band effective 
radii, which are plotted against Lj,oi estimated again from Lpjn (small triangles). Dashed line is the limiting luminosity given in Eq. 
(21). Two arrows at the bottom left corner indicate the effect of concentration parameter c on the resulting re- 



Actually, the limit on bolometric surface brightness we 
derive depends on the precise value of adopted concentration 

parameter c. In Fig. 14, we adopt c = 2.2 for all cases. As 
discussed in Section 4.4, it is possible that this concentration 
parameter could be large for UVSBGs. In Fig. 14, we show 
the correction factor for the intrinsic effective radius in the 
case of c=1.0 and 1.5 by arrows at the bottom left corner. 
Nevertheless, the derived limit on the surface brightness is 
robust for active starburst galaxies, since we confirmed in 
Section 4.4 that for ULIRGs and relatively luminous UVS- 
BGs the adopted value c=2.2 is reasonably consistent with 
observed effective radii. Note that even if less active star- 
bursts have the smaller value of concentration parameter, 
the bimodal feature of starburst activity shown in Fig. 14 
still remains. 



In summary, the SED diagnostics find that 1) UVSBGs 
have an intrinsic effective surface brightness LboirZ^ — 10^^ 
Lokpc^^, and 2) there is a limit on the size of the starburst 
region of around 0.3 kpc. Considering these results, we will 
propose a picture of bimodal starburst activity which con- 
sists of mild starburst activity in UVSBGs with a regulated 
surface intensity around L^oir^^ ~ 10^^ Lokpc"^ and an 
intensive mode of starbursts in ULIRGs confined in the re- 
gion of re ^ 0.3 kpc. In the next subsection, we will discuss 
the physical origin of this bimodality. 
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5.2 Origin of bimodal starbursts 

5.2.1 Disc instability m UVSBGs 

The regulating mechanisms for the mild starburst intensity 
have already been proposed by Meuror ct al. (1997). Fol- 
lowing Kennicutt (1989), they showed that the gas surface 
density is somewhat close to a critical value at which gas 
discs would bo unstable to self-gravity. According to Toomre 
(1964) and Quirk (1972), Kennicutt (1989) proposed the 
critical surface density for inducing star formation with the 
instability is 



Table 3. Velocity dispersions of starburst galaxies 



1^ 



(14) 



where ag is the velocity dispersion of the gas and k is 
the epicyclic frequency. The parameter a is a constant em- 
pirically determined as ~ 0.63. With angular frequency 
f2(r) = V(r)/r, the epicyclic frequency is given by 

1/2 



dr 



(15) 



While Kennicutt (1989) considered the instability 
around the fiat part of the rotation curve in disc galaxies, 
Meurer et al. (1997) applied it for the nuclear starbursts 
near the peak of the rotation curve in inner discs. We also 
apply the technique of Meurer et al. to interpret the bimodal 
trends of SFRs in Fig. 14. If we consider that the inner disc 
to be a rigid rotator with a angular frequency 0, the dy- 
namical time of the disc is given by 

With a constant Q, for the inner disc, the critical surface 
density is Ec = 2aa^l/{-KG). With the estimations for the 

dynamical time and the surface density, the maximum star 
formation intensity can be expected as 



9 X 10"^a(Tn^ Moyr"^kpc 



(17) 

^ and km s~^ kpc~^. 



where a and Q, are in units of km s 
respectively. 

With 'ip{t)/Li,oi from Eq. (11), the bolometric surface 
brightness limit can be introduced as 



LbolVe 



~ 3.1 X IQ^aon'^ Lokpc"^ 

= 3.7xl0^^f^U^^) 
Vo.aiy V20km s-v 



10'%m s kpc 



L0kpc 



(18) 



where we take t/to = 2 as a typical observable starburst age. 

Since most of disc galaxies have a central HI velocity dis- 
persion cTg ~ 15 — 20 km s~^ (e.g. Dickey, Hanson & Helou 
1990) and angular frequency Q, < 10^ km s~^ kpc~^ (e.g. 
Armus et al. 1990), this limiting bolometric surface bright- 
ness (~ 3.7 X 10^ Lokpc^ ) corresponds well to the limiting 
intensity of UVSBGs as seen in Fig. 14. Thus, starbursts in 
UVSBGs can be triggered by the self-gravitational instabil- 
ity in inner gas discs as already pointed out by Meurer et 
al. (1997). This instability seems to induce only mild star- 
burst activity, since the unstable region should have a scale 
of the disc thickness, and the whole gas in the disc cannot 
concentrate into the unstable region due to the centrifugal 



name 


<T [km/sec] 


Line"^ 


UBSBGs 


IC214 


124.6 


HI [21cm] 


JNG0114U 


62.3 


HI [21cm] 


NGC1569 


23.1 


HI [21cm] 


NGC1d14 


75.0 


CO 2.3//m 


NGC1d14 


107.7 


CO(l— >2) 


JN 0^4194 


41.7 


Mi [21cm] 


IN lor'o^iy^ 


1U4.U 




NGC4194 


73.2 


CO(1^2) 


NGC4385 


36.7 


HI [21cm] 


NGC5236 


78.2 


HI [21cm] 


NGC5253 


25.8 


HI [21cm] 


NGC6052 


127.8 


HI [21cm] 


NGC6090 


61.4 


HI [21cm] 


NGC6090 


50.0 


CO 2.3/im 


NGC6090 


55.2 


C0(1^2) 


NGC6217 


61.6 


HI [21cm] 


NGC7250 


52.0 


HI [21cm] 


NGC7552 


62.9 


HI [21cm] 


NGC7673 


51.1 


HI [21cm] 


NGC7714 


65.0 


HI [21cm] 


NGC7714 


106.6 


CO(l-+2) 


ULIRGs 


IRAS12112 


146.4 


CO(1^2) 


Mrk273 


71.8 


CO(1^2) 


Mrk273 


160.0 


CO 2.3/im 


IRAS 14348 


107.7 


CO(1^2) 


Arp220 


182.3 


CO(1^2) 


Arp220 


150.0 


CO 2.3/im 


IRAS22491 


85.6 


CO(l->2) 



Note: " Emission lines (HI 21cm; CO(l— »2)) and absorption line 
(CO 2.3/im) adopted to derive the velocity dispersion. See the 
caption of Fig. 15 for references to the observations. 



force from the disc's rotation. This means that the maxi- 
mum size of starburst regions is expected to be similar to 
the disc thickness ~ 0.3 kpc, and therefore consistent with 
the SED fitting results in which no staxbursts with re ^0.3 
kpc are found. 

The surface brightness of ULIRGs are an order of mag- 
nitude larger than that expected from the disc instability. 
Therefore, ULIRGs should be triggered by another mecha- 
nism which can induce a stronger concentration of gas. As 
indicated by Eq. (18), the surface brightness of starburst re- 
gions in UVSBGs are related to the physical properties of 
parent galaxies. On the other hand, the strong mass con- 
centration in the central region of ULIRGs should result 
in the strong self-gravity of starburst region, and therefore 
the starburst region in ULIRGs can be dynamically isolated 
from parent galaxies. In such a case, the surface brightness 
of ULIRGs should be controlled by the relative strength of 
the self-gravity and its feedback. 



5.2.2 Self-gravity vs. feedback in ULIRGs 

For a starburst with a typical duration of to, the kinetic 
energy of gas clouds with random velocity Vg can be written 
as l/2MgVg ~ Lkinto, where Lkin is the kinetic luminosity 
due to feedback. Assuming that Lkin is proportional to the 
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Figure 15. Comparison of velocity dispersions in starburst re- 
gions. We estimate the velocity dispersions from the SED fit- 
ting results; i.e, afit = {GM(< re)/2re)-'/^ where we assume 
M{< Te) = 2Af»{< Te). Squarcs and circles indicate ULIRGs and 
UVSBGs, respectively. For solid symbols, we adopt the velocity 
dispersions derived from the line profile of the CO(1^2) emission 
(Sanders, Scoville & Soifer 1991), while HI line width measure- 
ments arc used for small open circles, which are taken from LEDA 
database. To convert both CO and HI line widths to the veloc- 
ity dispersion, we assume that the line profile is Gaussian; i.e., 
CT=FWHM/2.35= W20/3.62, where W20 is the width at 20% of 
the peak. For large open circles and open squares, we adopt ve- 
locity dispersions derived from the 2.3 /im CO absorption feature 
(James et al. 1999; Shier & Fischer 1998). For open triangles, we 
plot terminal velocities, instead of velocity dispersions, which arc 
obtained for Mrk273 (450 km s'^), NGC1614 (330 km s"!) and 
NGC7552 (500 km s'^) by Heckman et al. (2000). 



bolometric luminosity with a constant efficiency fkin, i.e., 
Lkin ~ fkinLbol, we Can write the kinetic energy per a unit 

mass as 

1 2 ^kintdyn r ^bol /in\ 

By using the relation between Lboi and SFR tp in Eq. (11), 
the typical velocity can be written as 

Note that the effect of feedback is independent of the mass 
in the starburst regions and also in the host galaxies, even 
though a larger gas supply can induce a more active star- 
burst. On the other hand, the gravitational effect becomes 
more prominent in more massive systems. This means that 
too massive starbursts cannot be sustained by feedback 
against strong self-gravity. 

This limiting effect of feedback against self-gravity can 
be estimated by the comparison of Vg with the escape ve- 
locity Vesc,c {Ve3c,c — 2GM{< re)r~^), wMch indicates the 
depth of the gravitational potential well. A critical mass- 
radius relation is expected for the condition, Vesc,c = Vg. 
This relation can be rewritten to another luminosity-radius 
relation by using the mass to light ratio of starburst regions 
[Eq. (12)]. We can then introduce a limiting bolometric lu- 



minosity 




for ti, = 100 Myr, where we use M(< re) = 2M*(< re). The 
total mass within re is somewhat uncertain, since a realistic 
distribution of gas and the contribution from the dark mat- 
ter is important. Therefore, we simply assume A4{< re) = 
2M*(< re). At a typical starburst stage t/to = 2, we take a 
limiting bolometric luminosity Lboi.c — 1.7 x 10^^(re/lkpc) 
L0. The feedback effect cannot sustain the starburst re- 
gion against self-gravity when the bolometric luminosity ex- 
ceeds this limit. This limiting luminosity is indicated by the 
dashed line in Fig. 14, which crosses the line of limiting 
luminosity of the starbursts induced by the disc instabil- 
ity around 0.3 kpc. Fig. 14 shows that all the ULIRGs are 
above the dashed line where the feedback effect cannot sus- 
tain the starburst region against the self-gravity, while all 
the UVSBGs but one are below the dashed line. Therefore, 
the mode of starburst in ULIRGs can be different from that 
in UVSBGs as expected. 

The SED diagnostics for ULIRGs indicate a strong mass 
concentration in the starburst regions as shown in Fig. 13 b). 
Furthermore, the analysis of self-gravity and feedback allow 
us to specify that starburst regions in ULIRGs are too mas- 
sive to be sustained by the feedback effect. Such a system in 
ULIRGs should be dynamically unstable, which can be con- 
firmed by the observed kinematical signature. With the mass 
and effective radius derived from the SED fitting, we can es- 
timate the velocity dispersion as afu = {GM{< re)/2re)^^^ , 
which is independent of the direct kinematical estimation 
from line observations. If a starburst region is dynamically 
stable as a virialized system, the observed velocity dispersion 
should be close to this value. In Fig. 15, we compare afu 
with observed velocity dispersions. The observed velocity 
dispersions are derived from the emission lines of C0(1^2) 
and HI at 21cm, and the absorption line of CO at 2.3/jm. 
These velocity dispersions are tabulated in Table 3. The 
velocity dispersions by 2.3/im CO absorption are most suit- 
able to probe the gravitational potential in the starburst 
region, since the effect of bulk motion on the CO absorption 
feature is less important than that on gas emission lines. 
Comparing to the HI 21 cm lines, the emission line of CO 
molecules seems to be more suitable for our purpose, since 
these molecules are more concentrated in the starburst re- 
gion in general. However, we also rely on the measurements 
of HI 21cm line, since the observation of this line can be 
found in the majority of our sample galaxies, and HI line 
widths correlate well with the other kinematical tracers such 
as Ha and 0[II] which are strong in the starburst region 
(Kobulnicky & Gebhardt 2000). Note that cool gas is kine- 
matically decoupled with the outflowing gas which is indi- 
cated by triangles in Fig. 15. 

We can see in Fig. 15 that the sequence of the veloc- 
ity dispersion seems to be bimodal, again. The deviation 
of the observed velocity dispersion is most remarkable for 
ULIRGs. This means that, for ULIRGs, the whole starburst 
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region is dynamically unstable since the observed velocity 
dispersion is not enough to support it against the gravity 
due to large masses derived from the SED diagnostics. Such 
a dynamically unstable state for ULIRGs cannot be caused 
by the local disc instability as in UVSBGs, since it requires 
large mass concentration into the whole starburst region. 
Therefore, the dynamical disturbance of a galactic scale is 
necessary to activate ULIRGs. The merging of galaxies can 
be the most plausible trigger mechanism of ULIRGs. 



6 CONCLUSIONS 

We develop an evolutionary SED model of starburst regions 
in which the stellar and chemical evolution, including dust, 
arc taken into account consistently. Wo investigate the ba- 
sic properties of nearby starburst galaxies, such as the SFR, 
optical depth and apparent effective radius by using this 
evolutionary SED model. We studied the SED of various 
starburst galaxies with a SFR of 0.5 - 300 Mq yr~'^, clas- 
sified as UVSBGs or ULIRGs. Our SED model essentially 
reproduces all the SEDs of analysed starburst galaxies over 
the whole range of starburst activity. For starburst galax- 
ies having AGN, we find a clear excess of MIR emission 
when compared with the best-fitting SED model. The fit- 
ting results are confirmed by comparing the derived SFRs 
and optical depths with the emission line measurements. Ir- 
respective of the degree of reddening, our model can derive 
the SFR of starburst regions in a unified manner. The appar- 
ent effective radii derived from the SED fitting are found to 
be consistent with observations. Thus, our model can derive 
both the SFR and intrinsic effective radius from the SEDs. 

The variety of the starburst SEDs is caused mainly by 
the compactness of starburst regions, not from the starburst 
age. From the SED diagnostics, we find bimodal correlations 
between the SFR and the compactness/optical depth; more 
active starburst regions tend to be more compact and heav- 
ily obscured. 

The bimodal starbursts consist of mild state with a lim- 
iting intrinsic surface brightness Lt,oir~^ ~ lO'^^ Lokpc"^ 
and intense state with a characteristic scale of re — 0.3 kpc. 
The mild starburst can be triggered by the disc instability 
as proposed by Meurer et al. (1997). For ULIRGs, the sur- 
face brightness is an order of magnitude larger than that of 
UVSBGs. A simple analysis of the feedback allows us to de- 
rive a critical luminosity over which starburst regions cannot 
be sustained by the feedback against self-gravity. This lim- 
itation of the feedback causes the bimodality of starburst 
activity, in which the size of starburst region has a clear 
limit of Te ^ 0.3 kpc. According to this analysis, all of the 
ULIRG sample is identified as self-gravity dominated star- 
bursts. Furthermore, we compare the velocity dispersions in- 
ferred from the SED fitting results with observations to find 
that starburst regions in ULIRGs are dynamically unsta- 
ble. In order to produce the dynamically unstable starbursts 
with strong mass concentration, a violent trigger mechanism 
is required, rather than the disc instability. As suggested 
by various imaging observations, the most plausible trigger 
mechanism of ULIRGs would be galaxy merging which can 
cause the dynamical disturbance over a galactic scale. 
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